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An experiment using air-water mixtures to model steam- 
water flow in a reactor was devised to determine the turbulent 
mixing mass flow rates between two subchannels over a wide 
mimee OL flow conditions and therefrom predict turbulent 
ieernig OCCUrring in the reactor. An experimental apparatus 
was constructed for this experiment. 


Void fraction was chosen as a plausible correlating 
parameter for turbulent mixing and was measured together with 
axial pressure gradient for two subchannel geometries of 
interest. 


The subchannel geometries investigated were those found 
in a typical BWR fuel bundle. The flow regime in the sub- 
channels was detected by means of an electric conductivity 
probe and visual observations. Data obtained are compared 
with various correlations. 


Predictions of void fraction and pressure gradient with 
the Lockhart-~Martinelli correlation agreed closely with the 
Sa~ee) Comparisons were also made with Griffith's slug flow 
model, and with the Cravarolo-Hassid and Marchaterre-Hoglund 
wememiraction correlations. 
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CHAPTER I 


IMETOQUCE Ten 


The core of modern boiling water reactors is composed of 
a large number of subassemblies or fuel bundles. Each sub- 
meeemoly contains a matrix of vertical, Cylandrical, parallel 
fuel rods. The fuel bundles are generally enclosed ina 
Surrounding unheated wall. Figure:1l shows such a typical 
16-rod bundle arranged on a Square array. The main coolant 
flow is vertical, parallel to the axis of the subassembly. 
However, secondary radial flows between adjacent rods and 
between the outer rods and the wall also take place. 

The assumption is generally made that the phenomenon 
mewn as critical heat £flux (CHF), burnout, or departure from 
nucleate boiling (DNB) depends only upon the local conditions. 
The prediction, then, of the location and the magnitude of the 
CHF requires a complete and detailed knowledge of the main and 
secondary flow patterns. | 

Accurate estimation of the magnitude of the secondary 
flows, known as mixing, has been the major problem in the 
Peeaqiction of the critical heat flux condition. Considerable 
progress has been made in the area of mixing flows in single 
phase flow. The mechanisms of two-phase mixing, however, are 
still not well understood. 

It has been convenient to subdivide the main flow through 
the rod bundle or subassembly into flows between four adjacent 


rods or between a rod or rods and the subassembly wall. Such 
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artificially defined flow ducts are generally referred toes 
med channels or subchannels, 

The standard approach to a thermal-hydraulic analysis 
Sec) crOd bundle utilizes these subenannels 54 typical 
e@echannel breakdown of a rod bundle is ilidusteated in Figure 
1. Mean properties are assigned to each subchannel. Thus a 
feeieaxial velocity, enthalpy, void {or volumetric) fraction 
as assumed to exist over the entire cross section of a sub- 
femanel. This lumped PaLameter= approach, alenougn, fameumen 
Memo Strictly correct, serves as a first approximation to 
the actual velocity, enthalpy and concentration profiles ex- 
isting in a given subchannel. 

Mixing between parallel, ventilated (communicating) rod 
channels has been conveniently defined by Todreas and Rogers 
(2) as follows: 

Natural mixing 

ot) Ubbhalebeabeyae ibicereloveuars: = 

b) Diversion cross flow 
Forced mixing 

a) Flow scattering 

b) Flow sweeping 

iurbulent interchange is caused by the fluctuating 
velocities and pressures present in turbulent flow and has 
been broadly defined as a no net mass transfer phenomenon, in 
Ffeeerdance with Prandtl's mixing length theory. Such a 
definition has application to single phase flow but requires 


revision for two-phase flow. 


ae 
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Diversion cross flow is the net 2 nee flow resulting 
from a static pressure difference between parallel channels at 
a given level within the rod bundle. In general, it will be 
represented by a net mass transfer for single and two-phase 
Flow conditions. 

Forced mixing, by flow scattering and flow sweeping, 1s 
artificially induced. Rod spacers, for example, may produce 
a random diversion of flow (flow scattering). Devices such as 
helical rod wraps produce a preferential diversion of flow 
(flow sweeping). 

In the realm of natural mixing it has been assumed that 
the mechanisms of turbulent interchange and diversion cross 
flow are independent or at least may be analyzed as such as a 
first approximation. As noted by Todreas and Rogers (2), 
where large diversion cross flows exist, this probably ioran 
erroneous assumption, but nevertheless, it has been a 
reasonable method of attacking the mixing problem. 

The magnitude of the mixing flows has been investigated 
by Rowe and Angle (1), (8) and by Lahey, Shiralkar and 
Radcliffe (3) for the diabatic case. Rowe and Angle 
investigated turbulent interchange by using radioactive tracers. 
This method could not, however, give the net mixing vapor and 
liquid mass flow rates between two channels. It has been 
pew (e.g. (3)) that the vapor Bas san alftinicy vOmemigirace mae 
the less restricted regions of the central channels away from 


the unheated bundle walls. 


LS 








It would appear, then, that ewer oles emacs or tei oe 
the absence of a radial pressure gradient may be caused by at 
least two phenomena: 1) turbulent interchange, as in single- 
phase flow, and 2) a bubble drift to.less restricted flow 
areas. 

imespite Of this information, an) chvontne or omnia emacl 
mixing flow in the absence of a radial pressure gradient will 
@emeinue to be referred to as turbulent mixing.’ This as in 
concurrence with most literature. | 

mire mechanism Of diversion cross flow scems Lo Pewee iy 
well understood, and the magnitudes of these flows can be 
estimated by simply modeling the rod gap as an orifice subjected 
to a pressure difference. 

The mechanisms of turbulent mixing flows, on the other 
hand, are not at all well understood, and past Aetinataaee 
often been based on crude and somewhat arbitrary models. A 
typical method has estimated the turbulent mixing mass velocity, 
Gnixing’ as equal to some constant, 8, multiplied by the 
average mass velocity of adjacent subchannels. 


w!. (G. + 6) 
ae , i ik. 


Gnixing rc 2 
where wis is the mixing mass flow rate per unit length and c is 
the rod gap spacing. 

With this background, it can be seen that a more accurate 
determination of the magnitude and mechanisms aft turbulent 


mixing flows is required. 
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It was decided, then, to construceeon wa ppararica en 
could determine the magnitude of the turbulent mixing flows 
over a wide range of flow conditions. This apparatus, 
hopefully, could also give some insight into the mechanisms 
of two-phase turbulent mixing. Such an apparatus would have 
to exclude the possibility of diversion cross flow mixing 
and thus measure only the turbulent mixing component. This 
could be accomplished if pressure level and pressure drop 
could be equated at corresponding points in two adjacent 
feeenannels, thus removing the driving force for diversion 
cross flow. As a first approximation, such an experiment 
could be conducted under adiabatic conditions with air and 
water, modeling the steam-water flow of the reactor. 
Particular subchannels of interest are shown in Figure l. 


ee 


Subchannel geometry (b) was selected initially and is that 
geometry investigated in this thesis. 

It was furthermore postulated that void fraction might 
eetezgeswell as a correlating parameter for turbulent mixing. 
An indication of the influence of void fraction on mixing 
is the experimental evidence obtained by Rowe and Angle (8) 
and Lahey et al. (3), who showed mixing to be a function of 
quality, peaking somewhere in the slug flow regime. 

The purpose of this research work was the experimental 
determination of void fraction over a wide range of air and 
weeer mass flow rates for a typical rod bundle subchannel 


geometry. Simultaneously, pressure level and pressure drop 


were recorded ‘and flow regimes were identified. Thus a 


is e 








complete picture of the flow conditions was determined for 
a wide range of operating conditions. The data obtained 
have been compared with various correlations, and 


ppeitarities and differences are discucsce: 


1S 
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It 1s generally accepted that steam-water mixtures may, 
under certain conditions, be accurately modeled by air-water 
mixtures. In the present case this greatly simplifies the 
equipment necessary for oe study of two-phase turbulent 
interchange. 

As previously stated, turbulent interchange is defined as 
mre Mass transfer occurring in the absence of a static average 
pressure difference between adjacent subchannels at a given 
axial position. Hence, any experiment designed to determine 
this mixing flow would necessarily require equating pressure 
level and pressure drop at corresponding points in the sub-— 
channels. By establishing flows meeting these conditions with 
unequal void fractions in adjacent subchannels, one could 
observe the net mixing mass flow rates of air and water. 

The mixing flow rates could be obtained if some means 
were provided for Se Onereone the two subchannels upstream and 
downstream of a communication length, by simply measuring the 
flow rates of air and water in the two subchannels. This 
method is probably most accurate when the flows in the 
femacent subchannels have distinctly different qualities. 

The net mixing mass flow rates could also be determined 
meom tne difference in the inlet and exit alr and water mass 
flow rates,of each channel, with the channels continuously in 


communication. However, the maintenance of unequal void 


ley; 
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fractions and yet equal pressure and pressure gradients under 
these conditions is difficult if not impossible. Furthermore, 
in the first method, conditions of fully developed flow can 
ibe attained prior to allowing mixing to occur. This was 
thought to be a more realistic approach. 

The operating procedure then is to completely isolate the 
Subchannels, establish flows in each subchannel satisfying the 
conditions of equal pressure and pressure gradient, and measure 
the flow rates. Subsequently, a portion of the common boundary, 
the rod gad, is opened (that is, the subchannels are allowed 
to communicate over a variable length, the "window"). With 
the "window" open, the exit flow rates are measured again and © 
the net mixing flows of air and water are calculated. 


This net liguid mixing flow 1s then given by: 


- ' = sae ee ko se atl 
ee ee ee 
where 
wa = net liguid turbulent mixing mass flow rate per 
f 
unite Jength of channel (LB /HR ak) 
a a = turbulent mixing liquid mass flow rate per unit 
length from subchannel i to subchannel j 
(LB /HR ia) 
va a turbulent mixing liquid mass flow rate per unit 
f 
length from subchannel j to subchannel i 
(LB_/HR FT) 
os = length of communication (FT) 
ms ei ™. 6 = exit Mass TlowerareswOn led) i eoMm. tne 


ith and jth subchannels under isolated 


conditions (LB /HR) 


Wy 








a ° oe — 2 
lee a, il. = Same as m. - ; 
dTiegede ne ar ea 


are allowed to communicate 


m except that subchannel. 


A similar expression Gan be wreuttemmpen theme mets 


turbulent mixing mass flow rate. 


iy SG 


i — an 


m. =-(m.* -m. _) 


w' i ib oe _* - 
net,g Bee sre eae el i, See J1G Jeg 


Mmieenet total mixing mass flow rate:is oe = Eat + Uleeeee < 

The window length, L, will necessarily be chosen through 
experimentation. Hopefully, by measuring the mixing mass flow 
mtces with decreasing longths of communication, differential 
wemiges can be obtained at the limit by extrapolation to a zero 
window opening. 

BEwochannel geometry {5) Figure i, selected) for, 1m eral 
investigation, is representative of a typical wall subchannel 
Semee lett side of the rod gap. The Sight side of thewwed 
gap models approximately a central rod bundle subchannel, such 
as that existing between four rods. 

BeeeeocOgnized artificiality inthis expermenters seine 
effect See thic IMpOsStevon of a walleaeeene rod Se Ee isolate 
B@emenannels. Since the width of the rod gap (.168") is small 
[emeemparison to the total perimeter (2.071" and 2.301" for 
the left and right ehaniiom eect iTe ly). it may be argued 
meiae the effect on the overall flow conditions is small. 
However, the abrupt termination of ee wall at the beginning 
Sr the communication Window in the mixing apparatus will have 
Pome effect upon mixing. The eddies produced by the vorticity 


in the boundary layer undoubtedly will cause some additional 
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Mixing. It is argued, however, that with (sutticrvene larg 
window opening, the true turbulent mixing flows to be measured 
will dominate the artificially induced effects due to boundary 
ayer vorticity. 

A further effect of the presence of an isolating wall 
between the subchannels is the deviation of the void fraction 
feeeribution from that which would exist if no wall were 
present at the rod gap. It is expected, also because of the 
narrow rod gap, that this effect will be small. Nevertheless, 
if the two-phase Mixing length is found to be on the order of 
magnitude of the rod gap width, (1.e., if the average sub- 
channel properties are not representative of the local conditions 
where mixing occurs) concentration profiles in the communicating 
channels will probably have to be investigated. 


Sa le ee ie 


The net mixing mass flow rate per unit length for gas and 


- 


liquid combined can be postulated to be a function of the 


average void fraction, a, as follows: 


_— 


W W coe Wie 
nec gay 


cllas~os)oe Mis ¢ 


I 2 


J+ (a ,-a5)o,|u ia 


182 8 


In this formulation c is the rod gap spacing and |u. and 


i3,£ 
Mas g| are representative measures of the turbulent 
Fluctuating velocities of the liquid and gas respectively, 
eng. , RMS values or averages of the absolute values. 
The conservation equations and the turbulent mixing 


equation given above are developed in detail in Appendix A. 


21 





Since the net liquid and gas mixing flows are expected to 
ibe some function of the void £ractions ie che scoumun te eine 
channels, a detailed knowledge of void fraction is required 
over a wide range of water and air mass flow rates. In order 
to match pressure level and pressure drop at unequal void 
fractions, a knowledge of pressure drop as a iunctsonsone, oe 
fraction would also be desirable. Lockhart and Martinelli (5) 
have shown the void fraction and pressure gradient to be 
iMmetions of gas and liquid ee and hence, functions of 
pressure level. For this reason, the pressure level at the 
mecation of the communicating window is also of interest. 

Void fraction data for air-water and other two-phase 
mixtures in horizontal round tubes has been correlated by 
Lockhart and Martinelli (5) using a parameter, X, which is a 
function of flowing quality, Gensity and viscosity of the two 
phases. This same parameter, X, is also used to correlate the 


two-phase frictional pressure drop, dp/dl) to the single- 


els ENS 


phase gas or liquid pressure drop through the pressure drop 


faltiplier, on where 


m= N 


dl’ rR tp > ° 


mine subscripts T.P., 2, and g ee to two-phase, liquid 
Only, or gas only (respectively) filling the entire tube. 
Thus, for given liquid and gas mass flow rates, it is possible 
mempredict the void fraction, a, and the two-phase frictional 


meessure drop. 


Ze 





The elevation pressure drop 1S given by 
pg/G, = [pg (l-a)+ pial g/g, 


Hence, using the Lockhart-Martinelli correlation, the total 


pressure drop may be predicted. 


dp 


aoe = dp 
fe RO TAG, CP Seal 


FR TP * a1? ELEV TP 


ig 


) 


wore eraction correlations of other authors (12), (iS wae eco 
Meee in this work. . 

Another parameter of interest in two-phase flow is the 
ratio of average gas velocity to average liquid velocity, 


ffewn as the velocity ratio or slip ratio. 


As noted by Zuber and eae dase (9), in a flow channel, local 
Slip ratios exist as well as a local concentration (void 
fraction) profile. Knowledge of local void fraction may be 
mewpeni in the determination of a mixing length for turbulent 
mixing. Zuber and Findlay have defined a parameter, Cor 
which is a measure of the peaking of the concentration and 
velocity profiles. Average vo1d fraction data for a wide 
range of flow rates provides sufficient information for 


estimation of Co: 
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For round tubes, Zuber and Findlay assumed velocity and 


memcocntcrakioOn GLsSteri bu bvenseos 8 O lle ice 





R 

Jo 

OL, os at 
Cie = Ch R 

on WwW 


where j is the local volumetric flux density of the mixture 


(FT? /HR PT?) ana the subscripts ¢ and w refer to cond evens: 


the centerline and the wall respectively. 


The distribution parameter Co is defined as 





. 1/A f yUJGA aoe <a> ae 2 Gee tis 
O [1/A f,adA] + [1/A f,jdA] J. | oe m+n+2 


where the brackets, < >, denote cross sectional average values. 


i. nme Cvident that for adiabatic Elow (a. =a) 


E = ie ae op) ae 
O Tee aet 2 


Hor triangular concentration and VelLOCiiy. prOLriles 
t= n= 1), Cy reaches a maximum of 1.5. With values of m 


mmaem between 2 and 3 parabolic profiles are produced, similar 
to the velocity profiles for single-phase laminar flow. 


Bigher values of m and n (~7) produce profiles similar to 


Peebulent velocity profiles. As a flatter profile of 


concentration and velocity is reached, Cy approaches a value 


of unity. 
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Zuber and Findlay have also shown that the weighted mean 
gas velocity “a may be expressed in terms of o- and a weighted 


mean gas drift velocity V 


SJ 
= Sq) <aV > 
Vg = Ea eee Cy Vas 
where V = V -j) is the local drift velocity ef thavgae, 
OC) se 10 
ee |G £ 
<> = 


and 





= — <a 


Hence, C, can be determined from the slope of a i versus 
<j> plot, and the weighted mean drift velocity from the 
intercept with the ordinate axis. A typical plot in the 


e ~ <j> plane is given below. 





SLOPE =C, 





Hence it can be seen that plotting V, versus <j> over a 
wide range of flow conditions will yield, at least for round 
tubes, some information on the void and velocity profiles in 
meee rlow duct. Additional information on concentration 
profiles could be provided by a traversing void fraction probe, 
if necessary. 

The role of mixing length, &, in two phase turbulent 
mixing 1s of importance if one considers the net mixing mass 


How rate, w! asec LUMetL OMe Or war VOlLd teach 1Onmceu, 


net’ 
diffusivity, Ey! and a concentration gradient, oa 
ree 
4 28S = Snape ee ames 
Gaixing * “ce fy | Mys I eg hss | 
du 
= Ie. jaro le. | 
Lo, ( 4) 4! ay 
du 
" i 2 du 
_ 2 5 2 oi 
g 
= (p.-p_) 6 S% 
fg ady 
where 
gil. 
aS L dy 
and 
me | = appropriate measure of the turbulent mean fluc- 


Euabang velociiv es Of Inquldwandsqas, ) hece 


assumed to be equal (FT/HR) 
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— az 
So seine = net mixing mass velocity (LB /HR Pe) 
y = direction perpendicular to the main flow 
between adjacent rods 
u = local axial velocity 


The subchannel densities (94 7P5) nate refer to average 
representative values in the two aneenamele 

With knowledge of the mixing flow rates and void fraction 
profiles it may be possible to determine at least the order of 
Magnitude of 2, the two-phase mixing length. Although several 


investigators have proposed that G is independent of 


Mixing 
Gecmedap Spacing, Ss, and dependent rather on hydraulic diameter, 
Other investigators have proposed the opposite. Lahey and 
pemreauo (11) discuss this matter at length. The order o£ 
Magnitude of the mixing length, &, may shed some light on this 


—=— = 


SeomerOoversy . 

: While the above formulation seems intuitively correct, 
there are indications that it is probably true only under 
Sepeain conditions. One such condition might be the case of 
two identical subchannels in communication. In this case 
when the VOld| treactions are Caqual, sg = 0), one would expect 
zero net mixing. 

A broader view of this problem, however, suggests 
viewing the two subchannels in communication as a single 


Sheamnel progressing toward a condition of equilibrium 


Bemeentration and velocity distributions. 


ae 





When one considers a subchannel geometry such as geometry 


fal) in Figure 1, the, equilibrium Vold feacerenmat st cone om 
probably would not be one with ay ="0. “TAUS*suppoSii1 ones 


Supported by work discussed by Lahey and Schraub which shows 





that there 1S a distinct variation in equalibraum vyordstracrerom 
from the corner subchannel (geometry (a), Figure 1) and the 
centermost subchannels of the rod bundle. Such evidence 
Elpports the concept of a bubble drift mixing effect, 
previously mentioned in Chapter I. | 

ie ea and pee refer to average representative 


densities in the two subchannels under equilibrium conditions, 


then, using the same notation as before, one can show 


) |u;, 


Gnixing (P5-Ps og) ys to (P57 P sed | 


da 
= = -eq _ da 





Despite the fact that the influence of void fraction on 
turbulent mixing is not well understood, it would seem that a 
first step in any such investigation is the determination of 
average void fraction over a wide range of flow conditions 
for a given geometry. 

There are several experimental methods for measuring 
average void fraction. The simplest method is the use of 
quick closing isolation valves. If the valves can be shut 
simultaneously and SULEUC1Tently fast sto tiara weremegien Le EE 
sample of the two-phase mixture, the trapped liquid can be 


feeered to give accurately the average void fraction. The 
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application of this method to a non-circular duct necessitated 
the maintenance of the flow duct geometry throughout the 


valves in order to avoid perturbation of the flow. 
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CHAPTER III 


Description of Apparatus 


General 

The apparatus utilized to determine the required data 
was a basic air-water system with filters, vressure regulators 
and surge tanks for air and eae mixing devices, test 
sections, air-water separators and flow measuring devices. 
Other instrumentation provided were air and water temperature 
and pressure measurement, pressure drop measurement, and 
electric flow regime detection. The void fraction determination 
Was accomplished by guick closing valves and read directly 
through a transparent test section of lucite. Figure 3 is a 
schematic representation of the overall experimental apparatus. 
Letters in parentheses throughout this section represent_the 
various portions of the apparatus as indicated in Figure 3. 

Figure 4 is a photograph showing the major portion of the 
apparatus. On a support Framework ponciocine of Dexion 
strengthened by plywood was mounted the three foot entrance 
section (Q) and the five foot interchangeable test sections (S), 
Simulating reactor subchannels. The filters (B,F), pressure 
regulators (C;,G) and surge tanks (D,H) were located on the 
lower level with the mixing devices aN) VGCatedmca nee uilay 
under the lucite entrance section. 

Pe ene OxXte Of themvertical Cestascce ions a) aan Ooumene 
meper level of the laboratory, two identical air-water 


separators were mounted, one separator for each subchannel. 
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Photograph of Overall Apparatus 


Figure 4 
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Air flow measurement was accomplished by use of rota- 
meeers (X) at the exit of the aur Fromme use lam cols maa. 
air flow meters exhausted to the atmosphere. Because of 
meee wide range of aix £low to be measured, a bank of flow, 
meters with changeable floats was provided to meet the 
anticipated flows. 

tae Martinelli-Nelson correlation FOr Vold si rae-nuen 
versus quality was used to determine the steam-water void 
fraction corresponding to the maximum local quality existing 
in a boiling water reactor at 1000 psig. Using this same 
mema fraction for air-water, the Lockhart-Martinelli (5) and 
Wallis (4) correlations were used to predict anticipated air 
flow rates for a maximum water mass velocity of l x 10° I, Bm/ 
HRT“, PecduscsOtmene larger numboen "OG flowmeter > aac cuineacr 


—_ — 


air flow measurement was initially installed only at ene exit 
of the separators. Thermocouples (Y) were installed prior to 
the air flow rotameters. 

Water flow measurement, also by rotameters, was provided 
Beth prior to the mixing devices (K) and at the exit of the 
Peparators (2). 

From the beginning it was recognized that the inaccuracy 
of the rotameters (1.5% repeatability) might present problems 
in determination of the small mixing flows but would be 
adequate for void fraction determination. The advantages of 
a continuous rather than an intermittent device and the iuiltefial 


flow rates anticipated for the higher void fractions led to 


mae decision to utilize flow meters rather than an intermittently 
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operated receiving tank for flow measuring. Data using the 
mixing test section indicates however, that an intermittent 
flow measuring device will indeed be required, primarily 
because of fluctuations in the air flow at very low and very 


high flow rates. 


Separators 

The separators themselves used an inverted burette 
principle. The air-water mixture emanating from the test 
sections was introduced sp ehoromerseel Sayer eenl "Gah ~in -awlarcer 
memll tank (see Figure 5). The level in the outer spill tank 
was kept constant by the spillway, while the level in the 
inner can varied with air flow rate. The air-water mixture 
was separated by gravity and by inertial forces by use of 


—_ —-— 


Sarr les in the top portion of the inner "can". The water 
spilled into a spill. collection tank from which it was routed 
to an exit water flow meter. Slight variations in the water 


level in the inner can caused by the separation process were 


mioughnt to be the cause of the air flow fluctuations. 


Mixing Device and Air and Water Supplies 


iis Mixing device (N) 1S shown an detail in Pagune 6 
The air was introduced radially into the water stream by a 
Ziel pipe closed at the end and drilled with small holes 
along its length and circumference. The axis of the pipe was 
COincident with the axis of the water pipe and mixed flow pipe. 
The air supply chosen was the 200 psig supply in the 


Engineering Projects Laboratory. Since the mixing loop was 
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mee Only load on this line, the Supply weceeubyecc mare 
hmoenimum of fluctuations. If was found, hovever, jenatseecesuce 
of cycling of the compressor, it was necessary to bleed air 
from the surge tank for proper control by the air pressure 
fegulator at small flow rates. 

The water supply (street water pressure) was also subject 
to some small fluctuations despite use of the regulator and 
Surge tank. Data was taken largely at night to minimize these 
fluctuations. A water pump may be installed fOr Lucare 
Operation if higher pressures and less fluctuations are 
desired. 

Filtering was performed with a Cuno Micro-Kleen Filter 
fem~ealr and Fulflo Honeycomb Filter for water. 

Thermocouples (iL,4) were installed after the exits of 
the air and water throttle valves to monitor air and = 
temperatures prior to mixing. It was determined that these 
temperatures would best represent the air and water conditions 
in the test section. Only slight air temperature differences 


between the value before mixing and after separation verified 


this assumption. 


subchannel Sections 

mivemeesSt SECtioOn and preceding entrance section were 
machined in halves from lucite using a special milling cutter. 
Mee cntrance section and void fraction test section halves 
were joined with ethelene dichloride (plexiglass cement) to 


meeauce a strong joint cavable of withstanding the higher 


oi, 





pressures anticipated upon the closing of Ehe word freaeeren 
isolation valves. The mixing section halves were joined by 
bolts and sealed by O ring material. 

The entrance section and interchangeable test sections 
were joined together and to the inlet and exit blocks with 
rubber gaskets. 

The three foot entrance section was preceded by bypass 
memmections (0O,P) for use in the void fraction testing 
configuration. For least disturbance in the flow pattern 
the bypass connections were placed here rather than juste 
prior to the isolation valves. Thus a three foot long section 
identical to the test sections in geometry, was provided for 
establishment of fully developed flow prior to the test 
sections themselves. 


; An electric flow regime detector or flow vrobe (R) was 
installed in each subchannel. The TOCaEION OF Ene weer e cer 
was approximately four inches before the outlet of the entrance 
section. With the location of the flow probe in the entrance 
section, flow regime detection was possible for all test 
section installations. The flow probe itself was manufactured 
using a 24 gauge stainless steel spring cores Secon au £etlen 
Shield. A silver epoxy paste was painted on the lucite sub- 
channel wall for the wall electrode. Wire leads were imbedded 
in the subchannel wall underneath the silver epoxy to vrovide 


Emeurrent path. A detailed description of the electric flow 


Pw@e@oe installation is given in Appendix B. 
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Entrance Section 

The entrance section was constructed as indicated in 
Figure 7. A 1/32 inch brass strip was utilized as a subchannel 
separator, and this stripv was sealed in a machined groove in 
the entrance section wall. It would have been desirable to 
use a thinner strip in order to model more accurately the 
memedtions in the actual reactor rod bundle. However, two 
Semoicderations dictated 1/32 inch as the minimum strip thick- 
ness. Adequate strength under the possible pressure excursions 
from valve closing in the youd £fraciion test Section required 
a sufficiently strong strip. Secondly, the sealing of the 
bottom of the window in the mixing test section required the 
machining of a groove in the fixed (bottom) portion of the 
strip to mate with a knife edge in the upper, movable section 
memene Strip. It was this upper, movable portion of ene coe 
which allowed the adjustment of a variable sized window 
between the two subchannels for determination of turbulent 
mixing and for complete isolation of the two subchannels from 
each other when desired. 

Since the influence of hydraulic diameter on mixing was 
Meeorot concern, a "plug" was provided for the right hand 
Subchannel in the entrance and mixing test sections in the 
form of a brass flat bar machined to specifications and 
anchored by screws to the subchannel wall. With one set of 
brass bars in place, both the left and right subchannels were 
identical. With a second set of brass bars in place, the right 


Beechannel waS significantly larger than the left (see Figure 7); 
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Mixing Test Section 

The mixing test section was divided Anto upper vandelo wer 
Sections of three feet and two feet respectively, with a 
total length of five feet. The bottom of the window was 
Meecated 2 1/2 feet from the bottom of the mixing test section. 
The upper brass strip was movable to allow adjustment of the 
window opening and was withdrawn through a special packing 
gland in the exit block. A seal on the movable and fixed 
brass strips waS maintained by C clamp pressure on the exterior 


Seetche test section. Figure 8 shows the mixing test section. 


Void Fraction Test Sections 

A photograph of the void fraction apparatus is shown in 
Figure 9. The construction of the void Pneiion and entrance 
section from lucite necessitated a rugged mounting arrangement 
Merethne void fraction isolation valves. The void fraction 
test sections Sapeaene consisted of two 1 1/4 inch ball valves 
Semaeed On a 19 inch x 60 inch x 1/4 inch thick carbon steel 
plate. The plate in turn was bolted on the Dexion framework. 
Valve supports were provided by building a boxlike framework 
et 1/2 inch thick steel and 5/7 8 “anehwen els steel bored to 
accommodate the valve bodies. The valve bodies were machined 
at the inlet and outlet for brass sleeves which were soldered 
into the valve inlets and outlets. The sleeves served two 
Bmeposes: 1) to assist in fastening the valves to the boxlike 


framework and 2) to permit insertion into the valve ends of 
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Photograph of Void Fraction Apparatus 
Figure 9 





eylindrical lucite valve inseEets waren seoieay secnd= sic ea 

the same subchannel geometry as the test section. A similar 
cylindrical insert with convex ends was provided for each 
valve ball. Thus the subchannel geometry was maintained 
throughout the entrance section, isolation valves and void 
Mmeeetton section. The only interruption in the subchannel 
Seemetry were the ~.020 inch gaps at the valve balls to 
Peeerude hangup Of the ball inserts on the Buna N valve seata, 

Hiew approximately 43° 172° tneh SGI eon SecerOnebeeycem 
meee balls permitted vo1d fraction determination from 
meecoximately 3 - 97%. The 0 - 3 and 97 - 100% range of a 
was not capable of being read directly because of its being 
hidden by the valve bodies. 

The portion of the test section between the valves was 
feovrded with three joints sealed by O rings. Two of the 
joints were located at the juncture with the valve inserts 
and the third joint was located in the center of the section. 
Limitations on the traverse of the milling machine required 
the test section between the isolation valve inserts to be 
divided into two, approximately 20 inch long, sections. 

The void fraction test section apparatus consisted of a 
left subchannel and a right subchannel section, each capable 
of being installed while the main steel mounting plate was in 
place on the Dexion structure. Thus, there were interchangable 
valve inserts, valve ball inserts and void fractiGn ESOlatvem 


sections for both the left and right hand subchannels. A 
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photograph of these pieces 15 (shown anh gue Or ona een 
shows the valve ball insert in place. 

Void fraction divisions from 3 ~- 97.5% were scribed on 
the front only of the isolation sections by use of a milling 
machine. Visual error in reading void fractions due to 
refraction was found to be minimal sine the meniscus provided 
a gauge for determining when the line of sight to the air- 
water interface was above or below the line of sight 
perpendicular to the interface. 

Pressure taps for the pressure level and pressure drop 
determination were Gaste ted three feet apart in the isolation 
section of the void fraction test sections. Pressure trans- 
ducer and transducer readouts were utilized. This 
instrumentation permitted accurate determination of the 
Mie@eooure level and pressure drop. 

‘The valve closure device selected was an air SLL Sie ts 
operated mechanical linkage actuated by an electric solenoid 
Simultaneously with a similarly operated bypass valve. The 
main isolation valves were mechanically linked by a steel 
Operating rod and operated by a single air piston and cylinder. 
The bypass valve wasS operated by a similar linkage but with a 
Separate air piston because of the large distance from the 
test section to the bypass valve. The bypass valve fnd 
Operating mechanisn) was mounted on a separate Dexion and wood 
framework on the lower level. 

Fach air operated piston and cylinder was actuated by 


separate solenoids connected electrically in parallel. 
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Figure 10 
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Photographs of Valve Ball Insert in Valve and Valve Insert 
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Actuation of the solenoids closed the isolation valves and 
simultaneously opened the bypass valve. Careful adjustment 

@r air operating pressure to the separate cylinders gave a 
reliable simultaneous operation in 70 milliseconds. Preliminary 
Saleulations using the Wallis correlation (4) for air and 

water showed this to effectively isolate the section in ~302 

of the transit time for the gas phase at 80% void fraction 

and Bee ay = 1x10° lb/£t?hr. The enonel was effectively 
isolated in 1/8 turn of the isolation valves and completely 
isolated in 1/4 turn. Valve closing times were monitored on 


a dual beam storage oscilloscope using an electrical circuit 


with potentiometers mounted on the valve stems. 
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CHAP Ena 


Procedure 


#he right subchannel vo1d fraction) vest sscecerommerci tae 
with valves and inserts was mounted first above the right hand 
Miechannel of the entrance section. By so doing, the left 
subchannel of the ae Severs section was blocked off. The bypass 
valve was connected to the right subchannel bypass connection. 

The inlet water flow rate was set at predetermined 
Miigerticial liquid mass velocities matching the range of mass 
Petecities in a typical boiling water reactor. Six superficial 
miss velocities ranging from 0.5 x ho” Wom be Ome 10° LB_/HR Pr? 
imoren a 0.1 x Oo" increment) were investigated. 

For each liquid mass velocity, the air flow rate, 
monitored at the exit of the separators, was varied throughout 
miewtlow range Of the air flow meter bank. Several void 
fraction readings (at least three) were taken at each air flow 
rate and from this, an average void fraction was determined. 
ie was found that in the bubbly flow regime there was little 
Beateer in the data, while in the slug flow regime as many as 
five aie Siewieica Peomeoonoles were necessary to achieve a smooth 
void fraction curve. Under identical conditions individual 
void fraction readings varied by as much as 8% in the slug 
flow regime. 

For each flow condition air temperature was recorded both 
prior to mixing with water and after separation. Water 


memoerature was monitored prior to mixing only. 


3! 


Both pressure level and pressure drop in the test section 
were recorded for each point. A low-pass filter was provided 
petween the Ap transducer readout and the vacuum tube voit=-— 
Meter to dampen the effect of the pressure drapresciliations 
faethe siug flow regime. 

The electrical conductivity probe and visual observations 
were used to determine the flow regime boundaries. Polaroid 
photographs of the flow probe traces were taken for permanent 
mecord of ne Eker conditions. Several representative photo- 
graphs are shown in Appendix B, Figure B-2. These photographs 
enable the Fe EO of typical slug lengths. 

It was found that the design of the separators was 
Memetang the vOold fraction to ™about 65%. At these high air 
[mbewetates, increased pressure drop in the exit piping and exit 


oe +e 


air rotameters caused a buildup of pressure in the inner, 


‘iT 


Meeerted Gan Of Ehe separators. @ihis pressure excecded) cc 
hydrostatic pressure maintained by the water in the outer drum 
ead resuited in escape of the air from the bottom of the inner 
can. Redesign or modification of the separators will he 
necessary to reach higher void fractions. 

When the required data had been obtained for the right 
Memeo channel, the corresponding void fraction test section 
was removed and the operations were repeated for the left 
Subchannel. 

Aiteration of the valve closing speed over a range of 


100-160 milliseconds showed no change in the average void 


maction measured for typical conditions in the slug flow 





fmegime. Although anitial valve closure ie siwe romouming= 
order of 70 milliseconds, an average of 100 milliseconds was 
maintained throughout pemaaes taking period. 
The instrumentation wés calibrated as follows: 
aye alr rotameters:” initial calibration wi eae know ae ne 
volume flow rates using a gasometer 
b) water rotameters: initial calibration with known 
water flow rates using a weigh tank 
c) pressure level transducer: calibration prior to 
each data taking ruineusing : mercury manometer 
(normally daily) 
ae Ap transducer: electrical calibration prior to -cach 
GaGa teak time icunes im@usme Ilyana ai ly) 
e) air and water thermocouples: no calibration 
It is emphasized that throughout the range of flow 
conditions covered by this experiment, no attempt was made 
to maintain pressure level sometand ae loOcCaecilon ineere 


test section. 
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CHAPTER V 


Results 


The experimental results of this thesis are presented 
graphically in this chapter for the highest and lowest values 
Smesuperficial liquid mass velocity (1.0 x mee and 0. Same Mole 
LB_/HR ep?) , 

Data are presented in tabular noon for these superficial 
mass velocities and intermediate values in Appendix C. 

Void fraction data are plotted against Wat the alr mass 
Flow rate, which fe, Gi eos propomtiond | Vtesoqualiay tor tac 
low air mass flow rates investigated. Flow regime boundaries 
are indicated on these graphs. Figures 12-15 show these data 
Pwemeompare them to correlations by Griffith (12), Lockhart— 
Peamemmel)i (5) and Cravarolo and Hassid (13). The Griffith 
prediction applies to the slug flow regime only and is plotted 
for a round tube having the same hydraulic diameter. The 
Lockhart~Martinelli correlation was developed for slabie Zaye! 
water and many other liquids in round tubes and horizontal 
petabatic flow. Cravarolo and Hassid's data was obtained for 
argon and water or ethyl alcohol in several geometries: tubes, 
emnuba, and a seven rod bundle in a circular shroud. 

‘The pressure gradient was plotted against void fraction, 
since this information is considered to be of interest in 
establishing the proper conditions for the forthcoming mixing 
experiments. A prediction boy the Lockhart-Martinelli corre- 


Metion for total pressure gradient versus void fraction, is 
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included for comparison. Measured local conditions (i.e., 
pressure level and temperature at the level corresponding to 
the window location in the mixing test section) were used in 
the Lockhart-Martinelli correlation. It is possible, however, 
me predict void fraction and total pressure drop using this 
correlation without knowledge of the pressure at the window 
[eeation. This could be done by iterative calculations of 
pressure gradient and void BASS eae: Pressure gradient data 
is shown in graphical form in Figures 16-19. 

Average slip ratio is shown in Figures 20-23, Sie 
feenaterre-Hogiund SGorrelvativonm( 14), £07 Slip srs. showin 
these figures for comparison. This correlation was developed 
memporeaict velocity ratios in vertical channels of boiling 
water reactors. These researchers investigate several geo- 
metries in formulating their correlation (which makes no | 
distinction for differing geometry). Both adiabatic and 


diabatic cases were studied for steam and water in rectangular 


and circular channels. Air-water experiments were conducted 
in rectangular channels only. Rod bundle geometries were not 
investigated. 
<j 
Average gas velocity, Vo = ass, was plotted against 


everage volumetric flux density, <j>, in Pigures 24-27, in 
the form set forth by Zuber and Findlay (9). These plots 
permitted the determination of Co: the velocity and concen- 
jmeatton distribution parameter, and ak the weighted average 


gas. drift velocity. In all cases investigated, ee was found 


— 


35) 





to be very small, less than 1 FT/SEC. Values of on of about 


i 25 were obtained for a superficial liquid mass velocity, 


6 


Ger ef 0.5 x 10 in both subchannels, Shon tic wie heremacs 


mekocity, Ge = 1.0x Tain a Value Of JabouUrr 7 Um ter Cy was 
obtained. 
Intermediate superficial mass velocities for the left 


subchannel will be investigated in-subsequent work. 
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CHAPTER VI 


Discussion of Results 


The void fraction data obtained agree quite tavoraee, 
with the various predictions plotted in Figures 12-15. 
The Griffith Slug flew correlation (12) was developed 
using round tubes, annuli, and circular shroud bundle geometries. 


The void fraction under adiabatic conditions is given by: 


e) 
fe 30 +V, A 
a eee 1S) 
where 
ee 
Vy, = K,vgD -|- Ky x ) 


Dis defined differently for round tubes and other geometries, 


and Ky and Ky are also geometry dependent. 


me Cirst glance it would appear proper to utilize Criucivenee 


rod bundle curves for estimation of the constants Ky and Ke im 
Meeecorrelation. Griffith's work entailed nes use Of a 220. 
shroud diameter bundle with seven rods. Keeping the shroud 
diameter constant, he varied the rod diameter, thereby altering 
the hydraulic diameter. The bubble velocities Griffith ob- 
served were bundle averages, not the average subchannel 
velocities of interest in the subchannels tested for this 
Biesis. Hence the asymptotic values of Ky and Ky Lor Goune 
tubes were used in applying Griffith's correlation. When used 


in this manner, Griffith's correlation gave results which 


egeecd very favorably with the data. 
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For the right subchannel, Grrtt ven sicon se tae ve ite. 
round tubes produced void fractions somewhat less than those 
measured. In the left, more restricted subchannel, the 
Segeetith prediction yielded a higher yvordst vac tone ia ii 
observed. Intuitively, one would anticipate generally lower 
void fractions here than in round tubes because of retardation 
Of the liguid in sharp corners, etc. The fact that this same 
effect was not observed in the right subchannel is subject to 
speculation. 

Comparison of Griffith's slug flow correlation with 
Zuber and Finley's (9) formulation gives the following 
relationships: 


°g 


~ moe 
O gt Mgt AKy VgD+K., (= Fyn 


ae 
0-40 
(14K, ) (4A) i Gi) 


Or in Zuber and Findlay's terminology 


<n eee 


= Gg ee 
= (14+K, 7 <3>4+K, VgD 


<j 3 
Sf EES el 
Cre tae 
O GJ 
hence 
ak, and V_. = K, ¥gD 


The values of Cy obtained from the iE versus <j> plots (Figures 


mee) yield values of K. Of approximately 0220 "and™ 0.25. ae 


2 
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— 


intercept on the ordinate axis, ea! 1s difficulties tolspecims, 
exactly but seems to range from zero to about 0.5 FI/SEC:- 
Use of the asymptotic valves sion (eee hee ae oi Ky 


and Ky and equivalent pipe diameter for D in Griffith's 
correlation agrees very closely with ae Id Odema tl oMnacLe aimed 
meom the ce versus <j> plots, corroborating the assumption 
that Griffith's correlation for cirular pipes applies to the 
geometry investigated. 

Gricfith's correlation used in thisemanner prediees 
values for C. of 1.20 and ee of 0.34 and Onooert/ CC ote me 
Mert and right subchannels respectively. 

The values of Co of about 1.20 and 1.25 > obtained for seine 


two liquid superficial mass velocities of 1.0 x 10° and 


Heo xX 10° LB_W/ER ry? respectively are in consonance with a 
miatter profile expected from doubling the water mass velocities. 
ir Zuber and Findlay's definition of ca for round tubes is used 
and the concentration and velocity profiles are assumed to be 


identical, exponent values of m and n of 3(G_=0.5 x Hon 


LB_/HR 
Pr”) and 4(G,=1.0 Be es LB_/HR Pp?) are obtained for both 
Subchannels. 

The hockhart—-Martinel i. voudm@meact1 On wsGeate eran siden 
below the observed data in the region of intermediate void 
meactions. in the low and high void fraction regions the 
correlation seems to agree quite well with the data. This 
effect does not seem to be flow regime dependent at the 


tower end of the region of discrepancy, since it is present 


in the bubbly flow regimes for the right subchannel at 


BS 


Gy = 1.0 x mon LB_/HR rr? ana in the slug flow regime for the 
other flow conditions analyzed. As the annular flow regime is 
approached, the Lockhart-Martinelli prediction agrees more 
closely with the data and may even give higher void fractions. 
Slug-annular flow was observed in the, left subchannel at the 
higher gas flow rates investigated, but complete annular flow 
Was never reached. The fact that the Lockhart-Martinelli 
predictions fall below the observed VOld” EiaCea Omi ay eae ets 
than above as one might expect, is unexplained. 
mnie Cravarolo-Hassid correlation (13) %8or steam water 
4S included in the void fraction versus Ns curves for a 
comparison with their correlation and for steam-water and air- 
water mixtures. 
The pressure gradient versus void fraction curves (Figures 
16-19) agree quite favorably in shape with the curves predicted 
peepee Lockhart-Martinelli correlation, but are consistently 
below the predicted values. 
Scatter in pressure gradient data, particularly at 
higher void fractions in the slug flow regime can be attributed 
momiwo Sources of error: 
1) pressure drop fluctuations present in slug flow 
conditions, the effect of which could not be 
Ponpletely eliminated by the low pass filter, and 

2) fluctuations in liguid mass flow rate caused by poor 
control by the water pressure regulator. These, 
though small, have a significantly higher effect on 


PLGSoUGe Gradient tian om VYold tfactlommycuem 


IS 





fluctuations were more prevalent in the slug flow 
regime). 
Average slip ratio, plotted versus gas mass flow rate, 

a? in Figures 20-23, seems to be in rough agreement with that 
meeaicted by the Marchaterre-Hoglund correlation. For the 
right subchannel, the observed slip ratios fell generally 
below the predicted values. On the other hand, observed slip 
ratios in the left subchannel were consistently greater than 
predicted values. | 


Since the slip ratio is related to void fraction, 


Su. =e 
ia Ul (1~X) a 





one would anticipate generally higher slip ratios for lower 
void fractions with equal quality flows in different 
geometries. 

mine observed slip ratio, for the left, more restricted 
subchannel, then would probably be viet: than that predicted 
for a more open (circular or rectangular) geometry, as is the 
case. The fact that the data for the right subchannel falls 
below the predicted curve, ataleast for ally bute wtne haghe: 
void fraction region, seems reasonable, since it may be even 
less restrictive than the average of Marchaterre and Hoglund's 
geometries. 

seatter in the slip ratio data may be albtrelouted towne 
mact that the accumulation of errors in the void fraction 
measurement and flow measurement become magnified in this 


display of data. This is particularly true at lower void 


var 





fractions, where a small error in void fraction 2s multi ied 
beyond proportion in the slip ratios the £4ce tae ties ceartee. 
ie not uniform about a slip ratio Gf aporo tmatelin wine eee 
lower void fractions may be indicative of inadequate operation 
of the mixing device at a low gas mass CLOW saves. 

An analysis of several representative electrical 
conductivity probe traces in the slug flow regime yielded 
bubble lengths varying from 1/2 - 6 inches with as much as 
meimenes between bubbles. This data was obtaine@ for the 


6 Z 


mont subchannel at G. = 0.5 x 10 LB,/HR UP ales) Ce EG « 


£ 
iheebow Conditions such as theses a wide tscatlrer amr vod 
fraction data occurred. The test section length, when compared 
with bubble lengths, was probably too small to achieve an 
accurate statistical average in a single sample. 

A change in pressure at the level corresvonding to the 
window location in the mixing test section will obviously 
affect void fraction. An experiment to determine the magnitude 
of this variation was conducted by approximately doubling the 
gage pressure at the window location for flow conditions in 
the bubbly flow regime. A decrease in void fraction of 
ee eimate ly .04 was observed. 

It will be necessary to throttle the exit of either or 
both subchannels in the mixing experiments to produce equal 
pressure level on both sides of the window. It is expected 
that the effect on void fraction caused by such: throttling can 
be predicted using the Lockhart-Martinelli correlation. It 
will be necessary, however, to substantiate this by experimental 


results. 
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CHAPTER VII 


Conclusions 


An experimental apparatus was designed, constructed and 
yielded satisfactory data for void fractiom) COveree Wicca 
of flow Be tetiouce Higher void fraction regions beyond those 
covered (oa = .70) may be reached with modification or 
redesign of the air-water separators. 

The data obtained agreed quite favorably with various 
correlations, although some discrepancy was present. Notably, 
the void fraction data were generally somewhat higher than the 
Lockhart-Martinelli predictions in the intermediate void 
fraction range. The Griffith slug flow correlation agreed 
quite well with the data in most cases. 

Pressure gradient varied with void fraction in such a 
manner that it should be possible to match pressure level and 
pressure gradient under flow conditions of widely different 
void fractions in the mixing experiment. 

Plots in the We versus <j> plane using Zuber and Findlay's 
Permeation gave a rough indication of the peaking of the 
malocity and concentration Bronce for varying superficial 
liguid mass velocities. neeanied profiles should be provided 
by other means, if needed. 

The electrical conductivity probe was very useful in the 
identification of flow regimes and provided a means for 
estimating-slug an bubble lengths. The design of the wall 
Piecterode for a normally non-conducting channel wall proved 


HO be satisfactory. 


Le 





Flow conditions in the subchannel geometries investigated 
are probably not representative of rod bundle average conditions, 
but rather representative of individual subchannels. Indeed, 
the purpose of these experiments was the simulation of sub- 
channels. Preliminary indications suggest that the right 
subchannel, which models the flow region between four fuel 


rods, is representative of the actual geometry. 
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APPENDIX A 


Two-Phase Conservation. EQuaetOls moon Voom ae ahem 


A better understanding of the influence of void fraction 
On two-phase turbulent mixing may be achieved through examina- 
men of the relevant conservation equations. 

Several multichannel computer programs have been developed 
using these eguations to predict flow conditions and ultimately 
the burnout point for nuclear power reactors with ventilated 
parallel channels. The COBRA Code, (7), developed for single 
and two-phase flow, was Selected aS an Cxonpl cron someceiie 
attempts to solve the problem of mixing as applied to burnout 
prediction. 

The assumptions of the COBRA Code are as follows: 
a) One-dimensional, two-phase flow is described by two 
average velocities, a gas velocity, Vay sunkol ey IUanGie uc. 


velocity, V The. velocity 07. sl ema enor S=V,/V ei 


f° 
1S generally larger than one. 
b) A static pressure drop can exist between adjacent 
subchannels, resulting in diversion cross flow 
an In diversion cross flow, the cross flow enthalpy is 
that of the subchannel from which it originates 
qd) Turbulent mixing can transfer enthalpy but not mass 
Also implicit in the assumptions is the representation oya 


Subchannel variables, such as enthalpy, velocity, etc, by 


Single average values. 


sh 





Lahey and Schraub (11). suggested neeecaie ier emass 
transfer approach cannot describe two-phase turbulent mixing 
adequately, and they propose an equal volume interchange as 
a more reasonable assumption from mixing length theory. This 
assumption requires experimental verification. The small 
amount of mixing data collected to date using the mixing 
apparatus indicate that for parallel channels subjected to 
equal pressure gradient there may be a net turbulent mixing 
flow by volume and by weight. 
Another point raised by Lahey, sheen ang Radciitice mio) 
is that the cross flow enthalpy is not always the same as that 
of the donor subchannel. One explanation for this condition 
may be a different velocity for each of the two phases in the 
diversion cross flow. 
An example with two parallel and communicating Sr syenata els 
has been selected to illustrate the use of the conservation 
equations, taking into consideration the effect of these 
proposed modifications to the equations presented in reference 
COE 
The following additional assumptions hawe Jere ALIN Ss 
Memeenese equations: | 
Peek net turbulent mixing macs ploweratel can sexi seem 
oon ase flow 

b) The diversion cross flow liquid and vapor velocities 
may differ 

er Anwoaxial velocity, su, is defined in order to conserve 


thes monecncum OL Lne mi cure- 
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YU = V_ (1-X) a Me 








where v' is defined such as to pres 
bscure : 
oe 1 i, ke 
(l-a) 9, 106 


Srve thee momo ell sein mie Ge 


The following mass continuity, energy and momentum 


equations apply to the control volume shown below for the 


momentum equation. 
os 
eae 
tk 


Z 


F:dx 


u 


ee i a ee en ee a er — 
' 


CHANNEL 


Control Volume 


3) 5. 


Showing 


XL j+det, 
Dj +d; 


b, 
s(n wn, gui pal, u) 


CHANNEL 
J 


Momentum Pillisce> 





MASS CONTINUITY EQUATION 


dm. 
-~ = - striae\iaeine ies ee A ee 
dx Lj spe! ley 7 ak 
where either ee Vieng Or both are zero depending upon the 


Peemscence and direction Of diversion cross Elow. 


eee = CGe” = eGii(l-a. 5 eral oe 9 ju fi 


a5 ig Ly eee. mg Nias 


where ju. and |u. | are appropriate measures of the 


i,t! Lp, Gh | 
turbulent components of the liquid and gas velocities 
Pempendicular to the main axial flow (x direction), e.g., 


RMS values or averages of the absolute values. Changes) @r 
g 


maemvOid fraction, Qs within the length dx are neglected in 


Mae above equation. Similarly, 
Wo = ie = = ee Q. a, ) oly; eit CsPglYas gl 
but Juss gl Jus; gland ee eel 
hence 
Wi Mas ie : 
= Ore AG Oe 3) eG, gltlaj,-asdo [uss g 


The diversion cross flow terms, allowing for slip in 


the direction perpendicular to the main flow become 
"ye melo: .)0¢ ace ¢ 7 Be Bela 5 a 


for flow from subchannel i to subchannel ernie 
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Wao = Ci Ol oF bia + GQ. 
ae! L ( oe ale Pa" ya g! 


Fh 
= 

tote 

Fh 


mei. t Low from subchannel jce.. 
The appearance of the void fraction difference between 
degacent subchannels in the turbtillent mixing masa flow races 
is a strong indication that void {raction can be used as a 
correlating parameter for Mixing. If one considers the 
piieipated influence of void fraction on the fluctuating 
fearing velocities, Juss gl and Mie els 


void fraction is an important paiwumeter in mixing seems even 


the assumption that 


more conclusive. 
Such a correlation need not hecessarily be linear, since 
(as noted by Lahey and Schraub (11)) turbulent mixing seems to 
be strongly dependent on the two-phase flow regime, peaking 
somewhere in the slug flow regime. Judging from the BReencbann ces 


caused by the passage of bubbles of various sizes, this 


suggestion seems reasonable. 
ENERGY EQUATION 


d(mh,) 


dx ia Wt Lame Sa We pre igi ee Wa ihe! 


J Jto3 je gpa 

Where, as before, either one or hoth of the diversion cross 

flow terms may be zero, depending on the existence and direction 

eeemee Civersion cross flows. | 
Assuming a uniform void fraction dqicstmibucton.) Eno ms iln6 


energy transfer terms are 
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ee ee wee we 


' ae = 13 A 
Ws ne, cre a deghe [Usa ¢]+az0 h er 


ing g i 


= 
5 
I 


c{(l~a,)ogheluzys ¢ltajoghy| |] 


er oa 
15+g 
Bee tne diversion cross flow energy transfer eterms are 


eo) Oe Ue eee ae ee 


ee 
fy 2 


w..h. = c[(l-a. | eves eee + a.oih ue: 
ye U 5) ¢ eal bok j?gog ji,g! 


MOMENTUM EQUATION 


GPs 
is a ee ee 


= “owe Ue Ow ee WwW. . 'ULGX Sa, Ua emai) 
Dae 2 a eee 1 2 eae 


elgg Uae -|- fing Boks “+ m.U;) 


Here, as in the COBRA Code, £, and fi account Lor MON] 


D 
uniform velocity and void fraction (concentration) profiles 


in a subchannel. The static average density is defined as 
m= (lq. + O. 
ps (l-a.) 9, ae 


Expanding the axial momentum change term 


Senet = Me Als +D.QM samc (Go Ue Wea eerie Wee 
hae ae ee SY th ak 1 1) ee Hg) 


m. [v!dG.+G.dvi] + u.(w..-W..tw.. ow...) 
ne me Magee bana 8, Fea i? 7) 1 
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2 
u.dm.tA.G.dv, + UW. (Ws.-W. sue ee) 
1, a teal eens 


I 


A. Cav! + 2u. (w..-w. .tw. 
[eT a 5 lar |S Pais Fs 


Do 
Po £50963 


PRS are 


where £. is the friction factor obtained assuming that all the 


‘ flow is liquid and 
g dp/ax) pap 4 
mene two-phase frictional pressure drop multiplier. 


Bence the momentum equation ean be redweed to 


2 2 
ole ea iy 2 oh. Gah 7 
“a ~ 2gp,p., ‘S:? +o ax 
Gy i Vey c 
EriGtion pressume Goo aece leratlon pressure 
drop GWE Beoabe i) ne 
iF J ee ef (Wee ( boom ) ee ee ee te ae) 
eC g A. ieee i) jie Dy 
Ee ci 
gravity pressure momentum exchange terms 
ec OD ; 


- | ‘ = tw... | Se ; 
' Bale OES, a» ee (2u, E pus) ] 


ie is SS eEOnE that a EhoneughiesStudy son Ene  eEurouU lemme 
mixing phenomenon as it affects energy and momentum transfer 
must consider the concentration and velocity profiles. The 
lumped parameter approach may be justified bUbeete actuae 


determination of a valid mixing length will be dependent upon 
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not only the mean void fraction and mean velocity in each 


ehannel but also on Concentration and \yoloei esl. Eiaowtelomor 
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APPENDIX B 


Description of the Blectrical @Gonduer ia mere: 


The electrical conductivity probes wenme mounced near une 
discharge end of the entrance section for both left and paanlte (ole 
eeenannels, In this location the probes allowed devecyitomor 
flow regimes for each of the three test section configurations. 
Figure B~-l shows the flow probe installation. 

The most difficult problem encountered in the construction 
of the probe was to provide a conducting electrode at the wall, 
Since the wall material was lucite and non-conducting. Several 
methods were investigated, and a conducting silver epoxy paint 
was selected. Weds paste-like material (EPO-TEK 418H) was 
painted in a very thin coat on the lucite walls of the sub- 
Channels for a total subchannel length of six inches, centered 
at the flow probe. A portion of the conducting wall was formed 
by the brass strip separating the two subchannels. In the right 
subchannel an additional portion of the conducting wall was 
memmea by the "plug" (brass flat bar) inserted in order to 
feeeeee the hydraulic diameter of this subchannel. Outside 
connection with the wall electrode was accomplished by 
threading a small multistrand wire through a hole in the 
lucite wall and laying several strands in a shallow groove 
machined in the subchannel wall. The groove was then filled 
with the silver epoxy, covering the wire strands. This 
meecuccd a sufficiently smooth surface. Ihe epoxy was cured 


by careful use of a heat lamp, insuring that the temperature 
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ELECTRICAL CONDUCTIVITY PROBE 
FIGURE B-I 


was not sufficiently high to cause plastic flow or warping 
muehe Ihatenkwey, 

The probes themselves were cut from 24 gauge stainless 
steel svring wire, rounded at the probe tip and shielded with 
thin wall teflon spaghetti, except for an approximately 1/16" 
bare portion at the tip. The teflon shield was not wet by the 
liguid in the subchannel and consequently provided insulation 
between the probe tip and the wall See eeOke - The overall 
pe@ee Giameter, including the teflon.shield, was about 0.035, 
approximately 8% and 10% of the hydraulic diameters of the 
right and left subchannels respectively. While a smaller 
probe would have been more desirable insofar as Jeeqtesiglershcikole 
of the flow is concerned, it would not have had sufficient 
mixture. A packing gland (Conax MIC-062A) was screwed into 
an external mounting block and provided a seal for the probe 
as it passed through the subchannel wall. 

Mhe electrical circuit was identical in principle to 
that described in reference (15). A six volt battery was 
used to supply power. Shielded cable was used for connections 
to eliminate noise. 

Typical probe traces for various flow conditions are 
peeve in Figure B-2. These Polaroid photographs of oscillo- 
scope traces were taken for each of the data points recorded 
femmtcne right subchannel. 

In explanation of the photographs, signals near the top 


of the photographs indicate a current path through the water. 


De 





Signals at the bottom of the photographic wndareatealitemuie fon 


of the current path by the passage of an air bubble. 
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Photographs of Typical Electric Flow Probe Traces 
Figure B-2 
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RLGhe. Supehamiet Gy = 0m le LB /HR Py? 








W Paice o at Flow <a> No. of do ~ ip 
3 window _: Regime Samples dz G inlets feinles 
window 
roe 
ievey 
LB ' 
m LB pst S x 
aR psia z rp Ie F 
ik 
HES 6 25292 0.4470 S Oro q, ~ 417 0) 63-0 
Mimo 2/.23 0.1310 S - Oe 19 420 80.0 592 
feee> 26.20 0.1252 S +63) ies -403 $0.0 SoS 
oes «624.99 0.1195 S 7620 8 oe) 200 Soe 
meee 6 6L2321506—6lh6OL 1110 S oe 6 ENE! o0r0 ayo Nees, 
meo 21.41 O.1010 ~B-5 ae oell es 3 -406 0)..0 Does 
58) Zieod 0. 1045 B - 500 S Pech sy 20, 0 so ee 
wooo 20.72 0.0990 B moron 5 910 80.0 Do 
mer 20.42 0.0979 B OG 3 aero) 89.0 Se EneS 
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Right Subchannel G = Q.7 x Loe LB /HR pp? 








at | | : 
"g | Cae aie nae eal es “g inlet “£ inlet 
LOGE LOG. 
LB LB | 
a psia nak) tal wii i 
[m2 29.53 0.1430 S 0.651 6 Woo Wee 62,5 
Mems2 27.76 0.1345 S$ 0.611 6 ADE iL. el 
9.46 28.62 0.1375 S$ 0.612 6 430 81.5 Bil 
7.92 27.81 0.1345 S 0.616 13 OS BiG eal 
me, 28.82 0.1390 S 05586 7 “O33 Bile eal 
Mee «(27.93 021340 Bs O.5e4 9 .457 81.5 61 
me 24.41 0.1170 B O.522 5 4.416) 64.5 61 
[el 22.59 0.1082 RB 09.434 3 ea ces | wen 
1.42 20.79 0.1045 B 0.342 3 1374 81.5 61 
Mone 21-79 0.0995 B 0.272 3 4.377 81.5 i 
[eeser 19.80 0.0950 B 0.160 3 200 1,8 61 


We 19.50 0.0935 B Oo 7 3 - 403 oS: 61 
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Right Subchannel G = 0.8 x Le LB /HUR ep? 


"g seca ieste REgtite aoe poe cas = "9 inlet _ inlet 
Loc. Aare 
oad psia ni | psi op op 
HR ae : FT 
M50 30.44 .1475 S 0.628 6 0.493 74.5 63.5 
HOn50 29.58 .1405 S 0.609 7 0.463 84.0 59.5 
e494 428.98 .1375 S 0.580 7 0.465 84.0 59.5 
Bes «227.94 .1340 S @e Si hes 0.493 81.0 s01.5 
fers 27.22 .1307 S 0.534 7 0.476 81.0 59.5 
4.33 27.05 .1292 B-S 0.527 3 0.433 81.0 59.5 
pmeo. 25.52 .1223 B 0.506 3 0.410 81.0 59.5 
moS2 23.78 .1140 B 0.396 3 0.390 81.0 59.5 
moo) 22.38 ««.1065 B 0.304 3 0.386 81.0 Bo) oS 
mao) 21.87 .1042 B 0.248 3 0.381 81.0 5025 
"627 20.98 .1001 B 0.164 3 0.390 81.0 59.5 
.242 20.00 .0955 B 0.052 3 0.400 81.0 59.5 
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Right Subchannel G = 0.9 x 10° LB_/HR pp? 





ae 
"9 he ee neal > Ae a "9 inlet “£ inlet 
LOC. ele 
mae psia wn ese att oF 
HR eT ical 
Meas 30.69 0.1486 = Oa60G 7 Opes 4 5 6325 
fom 29.93 0.1429 5 Oe59'G 7 O24 757 oL ol 
eee )606hl 29.8606 600. L427 5 02554510 OF 5035 Ci: 61 
mee 29.82 0 1420 S 02566 6 0.483 83 edi 
fees «6228.28 = «20.1375 B 0.504 3 0.450 83 61 
fee 6226.93 0.1282 B 0.433 3 O77 417) G3 61 
imo 624.2 O22 so B O22 0 3 024203 > 553 61 
1.470 24.9 0.1189 B 0.317 3 0.403 83 Sil 
1.018 yee2oe OO. 1105 B AOS. S Os Gomes 61 
Weee 22.46 0.1071 B 0.146 S 02 5 6e hos 61 
Ome 21.32 0.1015 B 0.026 3 OES eo 61 
Bas hiearci talon at Ne = 50 LB_/HR 


Flow Regime Key 
Soe lug annular 
S) SUK: 
BP oepulboly-sluq transition 


B Bubbly 
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APPENDIX D 





List of Equipment 


ieessSsure Level Measurement 


Stratham Pressure Transducer Model UC2 


stratham Model [Clee UNLVersgal Transducer Readout 


Differential Pressure Measurument 


stratham 11639 PM 60 TC Pressure Differential Transducer 


Ramapo Strain Gage Adapler Model S6A 100B 


Hewlett Packard Model 412A pe Vacuum Tube Voltmeter 


Air and Water Flow Measurement 


Fisher~-Porter Flowrator;:; (various models) 


Peeeiraction Isolation Valvos 


Rockwood Figure 101 Bal} Valves 


(1 e/ 2es =e laUshaleg: valves; 57 4 ~ bypass Valve) 


Electric Conductivity Probe 


24 gauge stainless stee] SPEER Ire 
MIC-~062~-A Conax Benet mcs, 


MeO-TEK 418H - High Temporature, Electrically Conductive 


Epoxy, (lpoxy Technology, Inc.) 
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